Critical scaling of finite temperature QED 3 in anisotropic space-time 
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We investigate the scaling behavior of the critical temperature of anisotropic QED in 2+1 dimensions with 
respect to a variation of the number of fermions Nf. To this end we determine the order parameter of the 
chiral transition of the theory from a set of (truncated) Dyson-Schwinger equations for the fermion propagator 
formulated in a finite volume. We verify the validity of previously determined universal power laws for the 
scaling behavior of the critical temperature with Nf . We furthermore study the variation of the corresponding 
critical exponent with the degree of anisotropy and find considerable variations. 



(N 

o 

> 
O 

ov 
& 

>: 

ON. 

m ■ 

o: 



X 



I. INTRODUCTION 

The properties of asymptotically free quantum field theo- 
ries under variation of the number of fermion flavors are cur- 
rently investigated with great efforts. In particular the corre- 
sponding quantum and thermal phase transitions from a chi- 
rally symmetric to a conformal phase are attracting a lot of 
interest. Within QCD, these efforts are triggered by the de- 
sire to understand the origin and size of the conformal win- 
dow and the associated potential for applications for physics 
beyond the standard model, see Refs. d-d] and references 
therein. Another asymptotically free theory is QED in two 
space and one time dimensions. QED3 in its strong interaction 
version serves as a laboratory to explore theoretical concepts 
in the comparably simple Abelian framework. On the other 
hand, QED3 also has important applications as a potential ef- 
fective theory for strongly interacting fermionic systems like 
graphene JfJ and high temperature cuprate superconductors 
Jj,|5[]. These superconductors are characterized by d-wave 
symmetry of the gap function. The variation of the number 
of fermion flavors in these materials is tied to the number of 
two-dimensional conduction layers. Another particularly in- 
teresting feature of these materials is that they require a dif- 
ference between the Fermi velocity, vf, and a second veloc- 
ity v a related to the amplitude of the superconducting order 
parameter. This fermionic anisotropy can be as large as |7D 
A = up /va ~ 10. Furthermore, both velocities are also 
different from the speed of light, c s . Thus, in the effective 
QED3 theory all three Euclidean directions are different, lead- 
ing to an anisotropic formulation of the theory. Consequently, 
the critical number of fermion flavors Nf does depend on 
(c s ,vf , va ), see Refs. ifsl— TT^Tl for recent studies. The theo- 
retical study of this situation sheds light on the applicability 
of QED3 to the situation in the superconductor but also offers 
the possibility to explore anisotropy effects in general with po- 
tential applications to a description of the quark-gluon plasma 
in high temperature QCD. 

It is therefore also of great interest to introduce finite tem- 
perature into the anisotropic QED3 framework. Then, in addi- 
tion to the quantum phase transition at a zero temperature crit- 
ical Nf, one expects thermal transitions for fixed Nf when the 
temperature reaches an 7V/-dependent critical value. Particu- 
larly interesting is the region where the thermal transition line 
merges into the quantum transition point. From very general 
arguments, universal scaling laws for such a situation have 



been derived in Refs. 
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the critical value Nf of fermion flavors 



with a scale k 

at zero temperature and a (non-universal) parameter a 



The 



power law \Nf — iV/| _s ° with critical exponent 0o con- 
stitutes a universal correction to the well-known exponential 
Miransky scaling due to the running of the gauge coupling 
ITill . In principle there are two possible generic scenarios for 
the power law with |Oq| > 1 and |Oq| < 1, shown in Fig.Q] 
For I ©0 1 > 1> the phase transition line approaches the zero 
temperature axis with infinite slope, whereas for |O | < 1 the 
slope vanishes at the quantum critical point. In QCD, pre- 
dictions for 0n from the functional renormalization group ap- 
proach lfl3l [Ml are consistent with recent results from lattice 
simulations IU6I1 and result in |6o| < 1. 

In this work we investigate the corresponding behavior of 
QED3 under the presence of anisotropics with vf — va 7^ 
c s . As a tool to determine the relevant critical tempera- 
tures, we use the Dyson-Schwinger equations (DSEs) for the 
(anisotropic) fermion propagator of the theory. Due to the 
anisotropic setup, the three momentum directions have to be 
treated separately, which prevents the introduction of hyper- 
spherical coordinates in the Dyson-Schwinger equations. It 
turned out that a natural framework to deal with this situation 




FIG. 1. The generic thermal transition lines for the power law 

_ 1 

Wffi — Nf I with cases |O | < 1 and |6 | > 1. The quantum 
critical points for zero temperature are marked with filled circles. 
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obey the Clifford algebra {7^, 7^} = 2 (^„. We consider iV/ 
fermion flavors and the anisotropic fermionic velocities vp 
and ua which will be implemented via a factor gi tflv acting 
like a metric. The index i = 1,2 indicates the node in consid- 
eration. In the high-temperature superconducting (HTS) sys- 
tem, the different nodes correspond to different zeroes of the 
fermionic energy gap function placed on the Fermi surface of 
the HTS system 

The metric for the first node is given by 



FIG. 2. The diagrammatic representation of the Dyson-Schwinger 
equations for the fermion (a) and gauge boson propagator (b). Wig- 
gly lines denote photon propagators, straight lines fermion propaga- 
tors. A blob denotes a dressed propagator or vertex, whereas a dot 
stands for a bare fermion-photon vertex. 
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whereas the metric for the second node can be obtained 
from interchanging the fermionic velocities. It enters the La- 
grangian in the form of 



is the formulation of the DSEs in a finite volume, i.e. in a box 
with antiperiodic boundary conditions. This approach, origi- 
nally introduced in Refs. IU7l , [l8tl . has been used to study the 
finite volume behavior of isotropic QED3 as well as the be- 
havior of Nj at zero temperature in the anisotropic setup fl2tl . 
The results agreed qualitatively with corresponding ones in 
other frameworks iflflflTtl . thus underlining the feasibility of 
the approach. In this work we generalize the approach to in- 
clude finite temperature effects. As it turns out, we obtain 
values of I Go I that depend on the size of anisotropy. 

This work is organized as follows: In the next section 
we discuss the structure of the Dyson-Schwinger equations 
(DSEs) of anisotropic QED 3 at finite temperature as well as 
the truncation scheme we are using. We also summarize some 
issues concerning the formulation of the DSEs in a box. In 
section|III]we present our results and discuss the scaling prop- 
erties around the quantum critical point and comment on finite 
volume effects. We finish with a short summary and conclu- 



II. TECHNICAL DETAILS 
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The order parameter for chiral symmetry breaking is the chi- 
ral condensate which can be determined via the trace of the 
dressed fermion propagator 



r dx^ p r-im^ (T> ^ (0) 0)), 
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with inverse temperature j3 = l/T and spacial momentum 
p = {p\,p-z) T . The corresponding expression for the dressed 
photon propagator is given by 



SF,j(oj p ,p) = / dr 



D^{uj p ,p) 



dr 



(2tt) z 

(5) 



and the Matsubara frequencies in both equations read 



A. The Dyson-Schwinger equations in anisotropic QED 3 at 
finite temperatures 



uj p = 2nT 



n t + \ for fermions 



n t 



for bosons 
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The details of the formulation of anisotropic QED3 have 
been discussed in Refs. and summarized in Ref. lfl2Tl . 

Here we only state the most essential details, using the no- 
tation of Ref. 0]. The fermionic four component spinors 



with integer numbers rit- 

Diagrammatically, the Dyson-Schwinger equations are 
represented in Fig. |2(a)| and Fig. |2(b)| In Euclidean space- 
time at non-zero temperatures, they are explicitly given by 
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with the momentum fc M defined by the difference p^ — and the node index i = 1,2. The inverse bare fermion prop- 
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agator is given by S ](p) = Ztflvy/g^^Pp with p = 
(u) p ,p) and the inverse bare photon propagator is given by 
Dq 1 (p) = Z^p 2 P^ y {p) with transverse projector P T . The 
dressed fermion-boson vertex is denoted by T' 3 ( p^ , q M ). The 
renormalization constants Z\, Z2 and Z3 of vertex, fermion 
and gauge boson respectively are each defined by the ratio of 
renormalized to unrenormalized dressing function of the cor- 
responding one-particle irreducible Greens' function. Since 
QED3 is ultraviolet finite, all renormalization constants can 
be set to one without loss of generality. The Ward-identity 
Z\ = Z2 is then trivially satisfied. 

Since we consider anisotropic space-time, we have to take 
care of the individual dressing of the Dirac components of the 
fermions at the different nodes. In order to keep track of the 
general structure of the equations we use the shorthands lfl2ll : 



(9) 



and 



p^i := i (ujp,p ) p^ , (no summation convention !), 

(10) 

where A Mt i denotes the vectorial fermionic dressing function 
at node i. 

The anisotropic expressions for the dressed fermion and 
photon propagators then read 



s f\ ( w p:P ) = B i (u P ,p) +i ^fgT v iv Pn,i, 



(11) 



P 2 



IL^ (uJ P ,p), (12) 



where Bi denotes the scalar fermion dressing function at node 
i and (ui p ,p*) the vacuum polarization of the gauge boson 
field. 

Next we explain our approximation scheme for the DSEs of 
the fermion and photon propagator. One of the most employed 
schemes is the 1/Nf expansion, which led to some success in 
the early analysis of the chiral transition at large Nf and zero 
temperature IU91, l2(ll . In this approximation the vector part of 
the fermion propagator is just a constant, i.e. dressing effects 
in the fermion wave function are not taken into account. Later 
on it turned out that this approximation is too drastic and leads 
to the wrong behavior of the photon propagator in the chi- 
rally symmetric phase. In Ref. ll2lll it was shown that a self- 
consistent calculation of the coupled system of fermion and 
photon-DSEs beyond the 1/Nf approximation leads to non- 
trivial power laws of the photon and fermion wave functions 



in the symmetric phase with values of irrational exponents de- 
pending on the details of the approximation for the fermion- 
photon vertex. A nontrivial fermion wave function beyond 
1/Nf was also shown to be important to correctly determine 
the dependencies of the chiral transition on Nf in the pres- 
ence of non-zero anisotropics fl2ll . We therefore have to take 
into account these effects here as well. To this end we employ 
the same approximation scheme as detailed in Ref. |[l2tl . We 
use a form for the fermion-photon vertex that agrees with its 
Ward-Takahashi identity to leading order in the tensor struc- 
tures. This ansatz is given by 

lf(^,a,,,g)=V» ^^+^g) (13) 

without using summation convention and the fermionic mo- 
menta p, q at the vertex. In the gauge sector we are cur- 
rently not able to solve the photon DSE self-consistently for 
general anisotropies. We therefore approximate the photon 
self-energy by a model which is closely built along the self- 
consistent result of the isotropic case. In particular it takes 
into account the irrational power laws discovered in Ref. l2lll . 
The gauge boson vacuum polarization then reads 




^ _ 9j Pa 9jPS 



e 2 N f ( 1 pf 



l&v F v& \\/p^Pi+e- 2 p] +2k P; 2 
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with p 2 defined in Eq. (O and p = (u> p ,p). In Ref. 112 111 the 
value of n for the vertex truncation Eq. ( TBI at the critical 
number of fermion flavors Nj has been determined as n = 
0.0358. In principle, this number may depend on the fermion 
velocities in the anisotropic case. Nevertheless, due to our 
current inability to calculate this dependence we will keep k 
constant in all calculations. We also fix the gauge coupling 
e 2 = 1 independent of the value of Nf and temperature T. 
This choice is one of many possible ones that serve to compare 
theories with different numbers of fermion flavors Nf. 

Having specified our approximations (fT3l. (fT4l) . we now 
proceed with the fermion-DSE. In order to extract the scalar 
and vector dressing function of the fermion we take appropri- 
ate traces of the fermion-DSE and rearrange the equations to 
end up with 



Bi (oj p ,p)= Te 2 ^2 



Aa,i (U p ,p)=l 



■E 



d 2 q B l {uj q ,q)g^D llu {uj k ,k) 
(2n) 2 B l (co q ,q) 2 + ^) 2 ' 

d 2 q J 2 («M g^D^iuj^kyj -q^ g?»D Xu (u k ,k )1 A^^p) + A^{u qi q) 



Bi (u> q ,q) + ( q itll f 



(15) 



(16) 
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where again we have no summation for the external index 
fi. In the following, we will investigate the finite temperature 
equations in case of equal fermionic velocities vp = Va, re- 
ducing the fermionic vector dressing function to two instead 
of three different components and leaving us with only one 
sort of nodes. We therefore drop the nodal index i. 



B. The DSEs on a torus 

Considering anisotropy and additional finite temperatures, 
the computational demands increase significantly in compari- 
son to the isotropic vanishing temperature investigations, even 
despite the applied approximations for the gauge boson. We 
therefore continue our investigations by changing the basic 
manifold from three-dimensional Euclidean spacetime to a 
three-torus, which is well suited for the evaluation of the 
Cartesian sums emerging from the boundary conditions the 
fermion and boson fields have to obey. To be precise, only 
the time direction for the boson has to be periodic and for the 
fermion field antiperiodic. The spacelike components are cho- 
sen to show the same (anti-)periodicity as the time component 
of the according field for convenience. In order to guarantee a 
valid interpretation of the timelike direction as a temperature, 
we have to make sure that the aspect ratio of time- and space- 
like 'boxlength' is appropriate. We therefore distinguish the 
boxlength in coordinate space as L\ ~ 1/T in time direction 
and equal lengths L2 = L3 = Lx in spacelike directions. 
The momentum integrals are then also replaced by Matsubara 
sums, 

/<&<■■•> ^ k (17) 

The discrete fermionic momenta are then counted by q n = 
J2i=2 3(2"x/ Lx)(ni + l/2)ii, where represents a Carte- 
sian unit vector in Euclidean momentum space. For the pho- 
ton with periodic boundary conditions the momentum count- 
ing goes like q„ = I] l=2 3 (27r/Lx)(n i )ei- The resulting 
system of equations can be solved numerically along the lines 
described in Ref . l22tl . 

As already discussed in the course of previous investiga- 
tions U2.H2I, evaluations on a finite volume have a consider- 
able impact on the quantitative results for the critical number 
of fermion flavors, but hardly change the qualitative behavior, 
when considering anisotropics. Since the introduction of fi- 
nite temperature is conceptually very similar to the previously 
studied anisotropics, we expect similar findings for volume 
extrapolations and postpone detailed studies on this subject. 
Instead, we concentrate again on the variations of Nj, keep- 
ing in mind that its absolute number will change in infinite 
volume calculations. 
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Vf = V A 

FIG. 3. The critical number of fermion flavors, NJ : plotted as a func- 
tion of the fermionic velocities vf = va, and the temperature T. 



III. NUMERICAL RESULTS 

A. The critical Ny in the anisotropic case vf = va 

In this section we present our results on the finite tempera- 
ture phase diagram for the case of equal anisotropic velocities 
vp = ?;a with c s = 1. We obtain the phase diagram for 
the critical number of fermion flavors, Np from solving the 
set of Dyson-Schwinger equations dl5l ) and ( TToT ) on a torus 
of (N x = 39 2 ) x (N T = 23) points and a box length of 
e 2 Lx = 600. In time direction, the sum of Np = 23 Matsub- 
ara frequencies guarantee a reasonable aspect ratio between 
the spacial and temporal number of lattice points, necessary 
for finite temperature calculations. The order parameter for 
chiral symmetry breaking is the chiral condensate (tro de- 
notes the trace in Dirac-space), 

(M) = T W ^- 2 tr D S F {oj qi q) (18) 

or, equivalently, the value of the scalar fermion dressing func- 
tion B(ujq,q) at the lowest Matsubara frequency and two- 
momentum. For each given temperature and anisotropy we 
determined the critical number of fermion flavors Np 

The results of our calculations are shown in Fig. [3] The crit- 
ical number of fermion flavors has been determined numeri- 
cally at each dot in the diagram. Subsequently a Mathematica 
interpolation routine converted these values into areas of equal 
Nj in different colors separated by contour lines. The contour 
lines therefore have to be taken as approximate results from 
the interpolation procedure and only serve to guid e the eye. 
In agreement with our previous investigations IU2I1 we find a 
decreasing critical number of fermion flavors for increasing 
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FIG. 4. The critical temperature for the chiral transition as a function 
of the number of fermion flavors for the case vf — va = 0.8. 
Shown are our numerical results together with a fit of the power law 
l |20| > in the region Nf > 1.25, i.e. in the vicinity of the quantum 
transition. 



anisotropy at effective zero temperature. When we increase 
the temperature at fixed anisotropy we again find a decreasing 
Nf for all anisotropic velocities taken into account. It is note- 
worthy that for anisotropic velocities larger than the isotropic 
value, vp = ua > 1, this decrease happens at a compara- 
bly slow rate, while for vp = va < 1 the critical number of 
fermion flavors rapidly falls off with increasing temperatures. 



B. Scaling in the pseudo-conformal phase transition region 

There is a considerable number of works identifying the 
chiral phase transition in QED3 as of pseudo-conformal nature 
due to the dimensionful coupling constant a, see e.g. I2H 
I23I l24ll and references therein. The scaling behavior close to 
the phase transition relates QED3 with other strongly coupled 
theories such as many-flavor QCD that has been subject to a 
number of recent studies in the framework of renormalization 
group equations lfl3l - [l5l l25U28ll . These investigations have 
shown that the explicit form of the symmetry-breaking scale, 
ksB, is to zeroth order given by the expression 

k SB ~ ko8(Nf t0 - N f )\N c ffi - Nf\~^5 

xexp , (19) 

where ko denotes a scale-fixing constant and Nf is the criti- 
cal fermion number for chiral symmetry breaking at zero tem- 
perature. Furthermore, Oq represents the leading order term in 
the Taylor expansion of the critical exponent in | Nf — Nf \ 
and a a further parameter. Whereas the form of the scaling re- 
lation $1% has been derived from very general considerations 
and is therefore universal for strongly flavored asymptotically 
free gauge theories, the parameter a depends on the theory in 
question, see Ref. ifTill for details. A similar type of scaling 



law also applies to observables such as the critical tempera- 
ture. 

In a formulation of the theory on a finite volume it is ex- 
tremely difficult to assess the complete scaling law ( fl9l ). Ex- 
ponential Miransky scaling dominates ( fT9l only very close to 
the critical fermion number and therefore cannot be resolved 
within the temperature range available in such a setup. This 
is also the case for our formulation of the DSEs on a torufl 
In the following we therefore strive to identify the power law 
part of the scaling relation, which is expected to dominate fur- 
ther away from Nf until finally a value of \Nf — Nf \ is 
reached where scaling ceases altogether. In the following we 
seek to locate the region where the power law scaling 

T CP ~fc |JV; i0 -JV>r^ (20) 

is applicable^. We come back to this point below. 

As an example consider the case of vp = va =0.8 shown 
in Fig. [4] The power law d20b nicely fits the calculated points 
in the region 1.25 < Nf < 1.75. For Nf < 0.9 significant de- 
viations of the calculated critical temperatures from the scal- 
ing law occur indicating the finite extent of the scaling region. 
For Nf > 1.75, somewhere close to TVS we expect the expo- 
nential Miransky scaling finally to dominate over d20l ). This 
region can only be investigated in the (even more demanding) 
continuum formulation of the DSEs, which we will address in 
future work. 

The values for |0 O | and the corresponding critical fermion 
numbers iVS at zero temperature are given in Figs. [5] and 
[6] In general we observed that the size of the scaling region 
where the power law d20b can be found, depends on the value 
of 1 0o I . It is quite large for |0o| < 1, as can be seen in 
Fig. [4] but turned out to be much smaller for |0o| > 1 and 
consequently the extracted values for ©o in this case have to 
be treated with some caution. 

Fig. [5] shows that the anisotropic velocities provide a pa- 
rameter that qualitatively changes the scaling behavior close 
to the quantum critical point Nf . Shown are results for our 
standard torus of (iV| = 39 2 ) x (N T = 23) points and a 
spatial box length of e 2 Lx = 600, denoted by red dots. In 
order to assess the effects of finite volume we also determined 
three selected points (the crosses) on an even larger torus with 
(iV| = 59 2 ) x (N T = 23) and e 2 L x = 900, requiring 
a considerable increase in invested computational resources. 
Finite volume effects are clearly present and tend to amplify 
the results for the critical exponent. However, this quantita- 
tive impact is rather uniform with respect to changes in the 



The reason is a matter of scales. Close to N'j there is a large separation 
of the dynamical scale, i.e. the generated fermion mass, and the intrinsic 
scale given by the dimensonful coupling a of QED3. This separation of 
scales necessitates extremely large volumes which requires a continuum 
formulation of the DSEs [ 22] . Lattice calculations also have this problem. 
2 Strictly speaking, the relation )20t is only an upper bound for the chiral 
phase transition temperature since it is insensitive to local ordering phe- 
nomena due to Goldstone modes in the deep infrared 1 14, 29], which may 
yield corrections to j20) . However, since in our formulation at finite vol- 
ume and small lattice sizes such corrections cannot be resolved we stick to 
d20> as a good approximation. 
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FIG. 5. The critical exponent |Oo| plotted as a function of the 
fermionic velocities vf = da. Red dots indicate results obtained 
with a box length of Lx = 600/e 2 , whereas the crosses indicate 
results on the larger torus with Lx = 900/e 2 



of | Go | ~ 5 for the isotropic case indicates the presence of 
dominant exponential Miransky-scaling with only very small 
power law corrections. On the other hand, for anisotropies 
vp = va < 0.9 the power law corrections dominate. Unfor- 
tunately at present we cannot further quantify this issue since 
the exponential Miransky scaling re gion very close to is 
typically small (see Table 2 in Ref. IU5I1 for an estimate) and 
therefore not accessible in our formulation on a torus, as dis- 
cussed above. Nevertheless we may be able to qualitatively 
infer the size of the scaling region by going further away from 
Nj . Indeed, as discussed above, we find a larger scaling re- 
gion there for |0q| < 1, whereas for |0q| > 1 the scaling 
region is much smaller in agreement with the considerations 
above and in Ref. IT-ill . 

The zero-temperature critical number of fermion flavors 
(Fig. [6]) decreases for increasing anisotropies thus confirming 
the behavior already found in Ref. fl2Tl . Note that the zero- 
temperature fit of the critical number of fermion flavors is not 
identical to our previous results. The (small) deviations arise 
from the fact, that both calculations were performed on simi- 
lar, but not equally sized tori. 
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IV. SUMMARY AND CONCLUSIONS 

In this work we have determined the critical temperature for 
the chiral phase transition of anisotropic QED 3 as a function 
of the number of fermion flavors Nj for different anisotropies. 
To this end we solved a coupled set of truncated DSEs for the 
anisotropic fermion propagator together with model input for 
the photon self energy determined from the isotropic case at 
T = 0. We found a generic phase diagram with decreas- 
ing critical temperature with increasing Nf. The transition 
FIG. 6. The critical number of fermion flavors for zero temperature, line hits the N f -axis at the quantum critical point 2V£ . We 
Nf fi plotted as a function of the fermionic velocities v F = «a • investigated the critical scaling with Nf around the critical 

point and verified a universal scaling law determined on gen- 
eral grounds for stro ngly coupled asymptotically free gauge 
theories in Refs. lfl3ill4ll . A characteristic quantity that char- 
acterizes this scaling region is the critical exponent |0o|- We 
found that the anisotropic velocities provide a parameter that 
drives QED3 from a setting with |0n| smaller than one to a 
region with large values of the critical exponent |0 O | 3> 1. 
This result is interesting but preliminary, since the gauge sec- 
tor of the theory was not taken into account dynamically in 
our framework. Thus the full consequences of our findings 
can only be addressed when fully self-consistent solutions of 
the fermion- and photon-DSEs are available. Nevertheless, 
our present framework serves as a qualitative guide and con- 
stitutes a proof of principle for the feasibility of dealing with 
anisotropic systems in strongly interacting fermionic theories. 



anisotropy, a finding in a gree ment with similar estimates in 
the zero temperature case lfl2lPl We therefore conclude that 
our results serve very well to identify at least the qualitative 
behavior of |0o|(i>/). We find a maximum with a critical ex- 
ponent much larger than one in the isotropic case vp = va = 
1. For smaller anisotropies, 0.5 < vp = v a < 0.9, the 
critical exponent decreases rapidly to values smaller than one, 
I ©0 1 < 1, while for larger anisotropies 1.1 < vp = va < 1.5 
the decrease is slower. 

The size of the exponent 0o is a measure for the size of the 
power law corrections to exponential Miransky scaling close 
to Nj Q , as discussed in Ref. 11411 . These corrections are in- 
duced by the running of the gauge coupling determined by 
the photon polarization Eq. ( TBI) and found to be small for 
1 0o I 3> 1 and large for |0q| <C 1. Thus the large value 
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